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In  a  number  of  different  fields,  it  is  known  as  a  result 
of  a  priori  theoretical  analysis  that  a  sequence,  [un]i 
obtained  experimentally,  has  the  representation 


Even  if  the  dimension,  N,  is  known,  the  determination  of  the 
parameters  and  X^  can  be  a  matter  of  some  difficulty. 

If  N  is  not  known  initially,  the  problem  is  naturally  more 
comple- . 

The  purpose  of  this  note  is  to  outline  a  method  designed 
to  yield  the  value  of  N  without  the  determination  of  the  ck 
and  X^.  In  some  cases,  this  is  all  that  i3  desired. 
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OH  THE  SEPARATION  OP  EXPONENTIALS 
Richard  Bellman 


1.  Introduction 

In  a  number  of  different  fields,  it  is  known  as  a  result 
of  a  priori  theoretical  analysis  that  a  sequence, 
obtained  experimentally,  has  the  representation 

N 

(1)  un  -  ,1 


K)' 
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cke 
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Even  if  the  dimension,  N,  is  known,  the  determination  of  the 
parameters  c^  and  A^  can  be  a  matter  of  some  difficulty; 
cf .  the  discussions  in  [l]  and  [2^  .  If  N  is  not  known 
initially,  the  problem  is  naturally  more  complex. 

The  purpose  of  this  note  is  to  outline  a  method  designed 
to  yield  the  value  of  N  without  the  determination  of  the  c^ 
and  A^.  In  some  cases,  this  is  all  that  is  desired. 

2.  Recurrence  Relations 

We  start  from  the  frequently  used  fact  that  a  sequence 
such  as  that  given  above  can  be  generated  by  means  of  a  recur¬ 
rence  relation  of  the  form 


(1)  v ( t  +  N)  -  axv(t  +  N  —  1)  +  • • •  +  aNv(t),  t  -  0,1,... 

with  v(0) , v(l) , . . . , v(N  -  1)  prescribed.  Let  v(t)  be  a 
solution  of  (l),  and  consider  the  determinant 


(2) 


Ck(t) 


v(t)  v(t  +  l) 

v(t  +  1)  v(t  +  2) 


v(t  +  k) 
v(t  +  k  +  1) 


v(t  +  k) 


v  ( t  +  k  +  1) . . . 


v(t  +  2k) 


Assuming,  as  we  may,  that  in  (l.l)  c^  /  0,  i  •  1,2,.  ..,N 
and  that  A^  /  A^,  i  /  J,  it  is  easily  demonstrated  that 
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(?)  (a)  Ck(t)  -  0,  k  £  N,  t  -  0,1,..., 

(b)  C^ft)  -  C^O^V  0,  t  -  0,1,. ... 

(c)  C^(t)  /  0,  k  <  N  -  1,  for  large  t. 

These  facts  can  be  used  to  determine  N.  Whatever  Infor¬ 
mation  Is  available  concerning  the  reality  of  the  Xk  can  be 
used  to  shorten  the  work. 

3.  Prediction  Theory 

Once  N  has  been  determined,  we  can  focus  upon  the  para¬ 
meters  a^,  and  from  these  obtain  the  c^  and  X^.  Perhaps 
the  easiest  procedure  is  to  determine  the  values  of  the 
which  minimize 

T 

(1)  Q(a)  -  jo(ut+N  -  a!ut+N— 1  "“•**”  VV  ' 

for  some  choice  of  T.  Since  this  problem  is  a  discrete 
version  of  the  basic  problem  of  Wiener  prediction  theory,  (the 
Kolmogorov  version)  methods  developed  in  that  field  (for 
example,  those  of  Levinson  given  in  the  appendix  to  Wiener’s 
book)  can  be  used. 

4.  Discussion 

The  problem  is  actually  more  complex,  since  what  is 
measured  experimentally  is  not  un,  but  wn  “  un  +  rn'  where 
rR  is  a  randan  quantity.  It  may  well  be  that  in  this  case,  a 
reasonable  procedure  is  to  minimize  the  expected  value  of  the 
quadratic  form  appearing  in  (3.1),  and  to  choose  a  value  of  N 
which  minimizes  this. 
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